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ABSTRACT 

We use the recent CMB power spectrum measurement by the Maxima experiment (Hanany et al. 2000) 
to constrain the redshift of reionization z rc . This becomes possible by combining the CMB data with 
cosmological parameters from various independent measurements, including Big Bang Nucleosynthesis 
(BBN) and X-ray cluster data. Most notably, our results provide a robust lower bound on z lc . We 
find that z IC > 15 (8) at the 68% (95%) confidence level, unless the Hubble constant is larger than 75 
kms^Mpc -1 . 

Subject headings: Methods: statistical — cosmic microwave background — Cosmology:observations — 
early Universe — large-scale structure of Universe 



1. INTRODUCTION 

In CDM cosmologies the first baryonic objects form at 
redshifts as high as 50 (see, e.g., Haiman et al. 1996 and 
references therein). At some lower redshift z re the forma- 
tion of these and subsequent baryonic objects leads to the 
reionization of the universe by the combined effects of the 
fairly hard UV radiation from AGNs and the softer UV 
radiation from young, star-bursting galaxies and possibly 
also from a population of early, more uniformly distributed 
stars (the putative population III stars). Typical estimates 
of z Te for CDM models lie in the range z rc ~ 10 — 20 
(Haiman & Loeb 1998). 

Recently, various problems with CDM on galactic scales 
have lead to the proposal that dark matter is "warm" 
rather than "cold" (Moore et al. 1999; Hogan 1999; 
Sommer-Larsen & Dolgov 1999 and references therein). In 
WDM models the redshift of formation of the first bary- 
onic objects is smaller than for CDM and increasingly so 
with decreasing WDM particle mass. It obviously still pro- 
vides an upper limit to the redshift of reionization. Ob- 
servational lower limits on the redshift of reionization can 
therefore be used to place lower limits to the mass of the 
putative warm dark matter particles. 

Spectral observations of quasars at z < 5 show evidence 
for a relatively early epoch of reionization. The lack of 
complete absorption shortward of the Lya line (the so- 
called Gunn-Peterson effect; Gunn & Peterson 1965) im- 
plies z re > 4 — 5, possibly even z ro > 5.8 (Fan et al. 2000). 

Constraints on z TC can also be obtained from data on 
anisotropies of the cosmic microwave background (CMB) 
radiation (Griffiths et al. 1999). Compton scattering of 
CMB photons on free electrons leads to the suppression of 
the primary anisotropies on small and intermediate angu- 
lar scales and hence a characteristic change of the angular 
power-spectrum Cg. 

So far the CMB data have been used to place upper lim- 
its on the redshift of reionization z IC < 30 (Griffiths et al. 
1999). In this Letter we use the Maxima-1 and combined 
Boomerang and Maxima-1 data (Hanany et al. 2000; Jaffe 
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et al. 2000) to place astrophysically interesting upper and 
lower limits on z Te . The Letter is organized as follows: In 
Section 2, we motivate the choice of a number of cosmo- 
logical parameters prior to analyzing the CMB data. Sec- 
tion 3 briefly explains our method of analysis, the results 
of which are summarized in Section 4. Section 5 discusses 
those results and provides an outlook. 

2. COSMOLOGICAL PARAMETERS 

Apart from the redshift of reionization that we are in- 
terested in, a cosmological model is specified by a number 
of other parameters. Unfortunately, today's CMB mea- 
surements alone still provide rather broad constraints on 
them when they are jointly taken into account (an extreme 
example of the astonishing consequences can be found in 
Figure 1 of Tegmark et al. 2000). To obtain meaningful 
results, one has to use other observational data, or one's 
own prejudice, to fix some of the parameters beforehand. 

To begin with, we have to specify what sort of energy 
is present in the Universe. Usually, one takes into account 
five different kinds, namely vacuum energy Oa, space cur- 
vature Ok, baryonic matter Ob, hot dark matter Ohdm, 
and cold dark matter Ocdm, where the values denoted by 
O are the energy densities of the respective type relative 
to the critical density. On the basis of standard inflation- 
ary theory we assume that there is no space curvature, 
Ok = 0, and also that the initial power spectrum is scale 
invariant and no tensor modes are present (see, however, 
the discussion in Section 5). Furthermore, dark matter 
shall be cold or warm, Ohdm = 0, since the effect of mas- 
sive neutrinos will be negligible on the angular power spec- 
trum of the CMB (Dodelson et al. 1996) - but see below. 
Obviously, this means that only three types of energy re- 
main in the Universe. They are linked via the relation 
Oa + Ocdm + Ob = 1. 

The density in baryons Ob can be constrained rather 
nicely using Big Bang Nucleosynthesis (BBN) . Taking into 
account a broad range of recent measurements (Tytler 
et al. 2000), it seems reasonable to set the value of Vt^h 2 
to 0.019 ± 0.002, where h denotes the Hubble constant 
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measured in units of 100 kms _1 Mpc _1 . 

Furthermore, measurements on bright X-ray clusters 
(Ettori & Fabian 1999) provide a robust estimate of the 
baryonic fraction /b, the ratio of gas to total matter in 
clusters. Both observations (Bahcall ct al. 1995) and sim- 
ulations (White et al. 1993; Evrard 2000) suggest that this 
is more or less equal to the fraction of baryons versus all 
matter in the whole Universe. Translated into our cosmo- 
logical parameters, we have 

/b = o ^ = (°- 069 ± 0.012)/i~ 3/2 + 0.04. (1) 

lib + ' *CDM 

Finally, the Hubble constant is taken to be h = 0.65 ± 
0.10, as a plausible compromise between various recent es- 
timates (Perlmutter et al. 1999; Mould et al. 2000; Reese 
et al. 2000; Patel et al. 2000; Freedman 2000 and references 
therein). 

The possible contribution to the dark matter density 
from massive neutrinos with free-streaming mass scale 
much larger than the mass of large clusters of galaxies 
would not be accounted for by this method. Because of 
the small mass of such neutrinos (< leV) , the contribution 
to f^DM can be neglected. That Hhdm -C 1 is also a pos- 
sible interpretation of the Super-Kamiokande atmospheric 
neutrino data, as discussed by Tegmark et al. (2000). 

3. METHOD 

Recently, Hanany et al. (2000) published a measurement 
of the power spectrum of the CMB from a highly resolved 
patch of the sky. We use their published data throughout. 
We have also used the combined Boomerang and Max- 
ima data of Jaffe et al. (2000) for our analysis and have 
obtained results very similar to the ones presented here. 

The measured CMB data consists of ten estimates of 
the power spectrum Cg averaged over a certain I range. 
In order to compare the measured CMB data to a given 
cosmological model, we first calculate a numerical spec- 
trum with CMBFAST (Scljak & Zaldarriaga 1996). The 
model spectrum is then binned in the same way as the 
observed data, and the resulting averages compared to the 
data by calculating \ 2 ■ For measurements following Gaus- 
sian statistics, we have 

10 / . (data) a (model) \ 2 

f=s(^H' <2 > 

where A^ data - ) and /\( modo1 ) denote the average power in 
the ith bin for the Maxima data and the given model, re- 
spectively, and Gi is the measurement error of the data 
point. Actually, the probability distribution function of 
power spectrum measurements is slightly skewed in com- 
parison to a Gaussian, so we correct for this effect with 
the offset lognormal method by Knox et al. (1998); Bond 
et al. (2000). 

Assuming that the individual power measurements are 
independent, a confidence level of 68% corresponds to a 
X 2 of less than 11.5, while a x 2 of 18.3 or more allows to 
reject a model at the 95% confidence level. Therefore, we 
always show the contours of these \ 2 values below. To put 
it simply, we ask the question "Which models fit well?" 

Normally, however, one does not perform this simple \ 2 
analysis. Instead, errors are assigned to fitted parameters 
by maximum likelihood. Most importantly, this means 



that the error is determined in comparison to the best fit, 
and will in general be smaller. Again, one can summa- 
rize this approach in a simple question: "Which model fits 
best?" 

4. RESULTS 

When we apply all the constraints mentioned in Sec- 
tion 2, the only cosmological parameter that remains free 
is actually the redshift of reionization z rc . From the \ 2 
analysis, we obtain a confidence interval of 18 < z ro < 25 
at a level of 68% and a range of 16 < z rc < 29 at 95% 
confidence. Figure 1 compares the data themselves and 
the model that provides the best fit. This fit is clearly 
acceptable, having a \ 2 of 8.3 for ten degrees of freedom. 

As explained above, we use various observational con- 
straints on cosmological parameters. To illustrate the ef- 
fect of errors in their measurement on our results, we sub- 
sequently omit each of them in turn. This means that we 
have a second free parameter, in addition to the redshift 
of reionization. 

In order to illustrate the methodical point made in Sec- 
tion 3, in Figures 2, 3, and 4 we also show the confidence 
regions obtained from maximum likelihood analysis. It is 
interesting to see how much smaller they are compared 
to the regions obtained from the x 2 values. We interpret 
this as an indication that one should be cautious and use 
the larger regions, in order to not rule out perfectly valid 
models. 

Figure 2 and Figure 3 display the constraints on the 
redshift of reionization obtained when the baryon density 
n^h 2 and the baryon fraction /b, respectively, are allowed 
to vary. Of course, this variation somewhat broadens the 
allowed redshift range. However, the change is not very 
large, in either case the lower bounds on z ro go down to 17 
and 13 for the 68% and 95% confidence level, respectively. 

The situation changes somewhat when we let the Hub- 
ble constant vary around the adopted value of h = 0.65. 
Figure 4 shows the results. Despite the rather large un- 
certainty of h, we can still constrain z ro : Unless h > 0.75, 
which seems unlikely when combining various methods 
of Hubble parameter estimation, we can conclude that 
8 < z ro < 32 at the 2<r-level, and at the la-level we find 
that 15 < z IC < 27. 

5. DISCUSSION AND OUTLOOK 

We have used the Maxima measurements of the an- 
gular power spectrum of the CMB to constrain the red- 
shift of reionization z lc . We found that astrophysically 
interesting bounds can be found when one constrains the 
allowed range of cosmological parameters with other ob- 
servations, including Big Bang Nucleosynthesis and X-ray 
cluster data. Even when the possible uncertainties of these 
parameters were taken into account, we still obtained a ro- 
bust estimate 8 < z IC < 32 at the 95% confidence level, 
unless h > 0.75. The implications of this result for struc- 
ture formation theory and the mass of the putative warm 
dark matter particle will be discussed in a forthcoming 
paper. 

One caveat is our omission of the possible contribution 
from tensor modes: tensor modes (gravity waves) con- 
tribute extra power on super-degree scales (£ < 100), so 
that when the total power (scalar and tensor modes) is nor- 
malized at low I, the scalar modes become lower and thus 
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the high-£ power looks suppressed. That effect is similar to 
that of reionization, but fortunately this degeneracy can be 
broken with future observations of the polarization angular 
power spectrum. For an electron scattering optical depth 
of r = 0.1 - corresponding to z ro « 13 for standard cos- 
mological parameters - the polarization signal should be 
easily detectable by the MAP and Planck missions (Ben- 
nett et al. 1995; Bersanclli et al. 1996), and possibly also 



the forthcoming new Maxima and Boomerang missions. 
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Fig. 1. — The best fit (red line) to the Maxima data (black crosses). The length of the crosses' arms in the x- and {/-direction indicate the 
width of the I range and the measurement errors of the power, respectively. 
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Fig. 2. — Likelihood contours for varying baryon density fl^h 2 . The Hubble constant is fixed at h = 0.65 and the constraint (1) on the 
baryonic fraction /j, = is in effect. The solid blue lines indicate the one- and two-sigma error regions determined from the x 2 values. The 
dotted red lines show the same for the maximum likelihood analysis. The shaded bar indicates the observational la range of O^h 2 obtained 
using BBN. 
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Fig. 3. — Likelihood contours for varying baryonic fraction f b . The Hubble constant is h = 0.65 again, but now Vl b h 2 is fixed at 0.0f9. 
Line styles are the same as in the previous Figure. 
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Fig. 4. — Likelihood contours for varying Hubble constant. Both /j, and Q^h 2 are now constrained as explained in Section 2. Again, solid 
blue and dotted red lines indicate the error regions for the x 2 analysis and the maximum likelihood method, respectively. 



